Relativistic corrections to transition frequencies of Fe I and search for variation of the 

fine structure constant 



O 

o 

(N 

> 
O 

00 
(N 



43 . 

Ph. 
i , 



: 

i/5 ■ 
o ■ 

• i-H ' 

1/3 ■ 
>Y 

43 : 

Ph. 



> 
00 
(N 

^1- 



O 



V. A. Dzuba and V. V. Flambaum 

School of Physics, University of New South Wales, Sydney 2052, Australia 
(Dated: February 2, 2008) 

Relativistic energy shifts of the low energy levels of Fe have been calculated using the Dirac- 
Hartree-Fock and configuration interaction techniques. The results are to be used in the search for 
the space-time variation of the fine structure constant in quasar absorption spectra. The values of 
the shifts are the largest among those used in the analysis so far. This makes Fe a good candidate 
for the inclusion into the analysis. 

PACS numbers: PACS: 31.30.Jv, 06.20.Kr, 95.30.Dr 



I. INTRODUCTION 

Theories unifying gravity with other interactions sug- 
gest a possibility of temporal and spatial variations of 
the fundamental constants of nature; reviews of these 
theories and results of measurement can be found in 
Refs. [H, Strong evidence that the fine-structure 

constant might be smaller about ten billion years ago 
was found in the analysis of quasar absorption spec- 
tra [1, H, H, @, H Hi- This result was obtained using 
the data from the Keck telescope in Hawaii. However, 
an analysis of the data from the VLT telescope in Chile, 
performed by different groups [3, [13] gave a null result. 
There is an outgoing debate in the literature about pos- 
sible reasons for the disagreement (see. e.g. [TTI. fl2l]). 

All these results were obtained with the use of the 
so called many-multiplet method suggested in Ref. [HI]. 
This method requires calculation of relativistic correc- 
tions to frequencies of atomic transitions to reveal their 
dependence on the fine-structure constant. Atomic cal- 
culations for a large number of transitions in many 
atoms and ions of astrophysical interest were reported 
in Refs. HE HE EE El El HE M HI HI- In present 
paper we perform similar calculations for neutral iron 
which was never considered before. Until very recent we 
were unaware about any lines of neutral iron observed 
in th e q uasar absorption spectra. It was Prof. P. Mo- 
laro [24[ who brought to our attention the existence of 
such data and the intention of his group to use them in 
the analysis. 

Using Fe I in the search for variation of the fine struc- 
ture constant has several advantages. First, the values 
of the relativistic energy shifts are large due to relatively 
larger nuclear charge (Z=26) and strong configuration 
mixing. Second, these values vary strongly from state to 
state which makes it hard to mimic the effect of varying 
fine structure constant by any systematics. Finally, iron 
is very abundant element in the universe. Great part of 
the previous analysis was performed using the data from 
Fe II. 

Calculations for Fe I are difficult due to large num- 
ber of valence electrons and strong configuration mix- 
ing. Its ground state configuration for outermost elec- 



trons is 3d 6 As 2 which is strongly mixed with the 3d 7 4s 
configuration. There is strong configuration mixing for 
the 3d 6 4s4p and 3d 7 Ap odd-parity configurations for the 
excited states. With eight valence electrons and strong 
configuration mixing full-scale accurate ab initio calcula- 
tions for Fe I would require enormous computer power. 
We have chosen a different approach. Below we report 
a simple method which is specially designed for systems 
with strong mixing of several distinct configurations. It 
combines ab initio Hartree-Fock and configuration in- 
teraction (CI) techniques with some semi-empirical fit- 
ting and gives very reasonable results at very low cost in 
terms of computer power. The approach is similar to the 
well-known multi-configuration relativistic Hartree-Fock 
method (see, e.g. (25[) and can probably be considered as 
a simple version of it. The accuracy for the energy lev- 
els of Fe I is within few per cent of experimental values 
while estimated accuracy for the relativistic energy shifts 
is on the level of 20 to 30%. Due to strong configuration 
mixing the results are sensitive to the distances between 
energy levels. Therefore, special care has been taken to 
reproduce experimental positions of the energy levels. 



II. METHOD 

It is convenient to present the dependence of atomic 
frequencies on the fine-structure constant a in the vicin- 
ity of its physical value cvq in the form 



ui(x) = ui n + qx, 



(1) 



where u>q is the laboratory value of the frequency and 
x = (a /ceo) 2 — 1 , q is the coefficient which is to be found 
from atomic calculations. Note that 



duj 
^ dx 

To calculate this derivative numerically we use 
w(+S) - u)(-6) 



28 



(2) 



(3) 



Here 6 must be small to exclude non-linear in a 2 terms. 
In the present calculations we use 8 = 0.05, which leads 
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TABLE I: Even and odd configurations of Fe and effective 
core polarizability a p (a.u.) used in the calculations. 



Set 


Parity 


Configuration 


a p 


1 


Even 




0.4 


2 


Even 


id 7 As 


0.4192 


3 


Even 


3d 6 4p 2 


0.4 


4 


Even 


3d 8 


0.465 


5 


Odd 


3cf 4s4p 


0.39 


6 


Odd 


3d 7 4p 


0.412 


7 


Odd 


3d 5 4s 2 4p 


0.409 



to 



10 (w(+Q.05) - w(-0.05)) . 



(4) 



To calculate the coefficients q using ^ , a must be varied 
in the computer code. Therefore, it is convenient to use 
a form of the single electron wave function in which the 
dependence on a is explicitly shown (we use atomic units 
in which e = K = l,a = I /c) 



III. 

iag v (r)fl(n)ji ri 



(5) 



where n is the principal quantum number and an index 
v replaces the three- number set n,j,l. This leads to a 
form of radial equation for single-electron orbitals which 
also explicitly depends on a: 



dfv K v j, , % 

~r + —fv{r) 

ar , r 



dg v 
dr 



2 + a 2 (e v -V HF ) 9v(r)=0 



-fv(r) + (e v -V HF ).f v (r)=0, 



(6) 



here K = {-l) l+3+1/2 {j + 1/2), and V HF is th e Har tree- 
Fock potential. Equation (|6|) with a — a^y/S + 1 and 
different Hartrec-Fock potential Vhf for different config- 
urations is used to construct single-electron orbitals. 

Tabic U lists configurations considered in present work. 
First four are even configurations and other three are odd 
configurations. We perform self-consistent Hartree-Fock 
calculations for each configuration separately. This al- 
lows to account for the fact that single-electron states 
depend on the configurations. For example, the 3d state 
in the 3d 6 4s 2 configuration is not the same as the 3d 
state in the id 7 As configuration. In principle, it is possi- 
ble to account for these differences in the CI calculations. 
One would need to have a complete set of single-electron 
states and construct many-electron basis states by re- 
distributing valence electrons over these single-electron 
basis states. Then actual many-electron states are found 
by diagonalization of matrix of the effective CI Hamilto- 
nian. This approach works very well for the case of two 
or three valence electrons (see, e.g. [H, H3, HH). How- 
ever, for eight valence electrons it would lead to a ma- 
trix of enormous size making it practically impossible to 
saturate the basis while using limited computer power. 



The results with unsaturated basis are very unstable and 
strongly depend on where the basis is truncated. There- 
fore, we prefer to account for the differences in the con- 
figurations on the Hartree-Fock rather than CI stage of 
the calculations. 

The self-consistent Hartree-Fock procedure is done for 
every configuration listed in Table U separately. Then va- 
lence states found in the Hartree-Fock calculations are 
used as basis states for the CI calculations. It is im- 
portant for the CI method that atomic core (Is 2 . . . 3p 6 ) 
remains the same for all configurations. We use the 
core which corresponds to the ground state configura- 
tion. Change in the core due to change of the valence 
state is small and can be neglected. This is because core 
states are not sensitive to the potential from the electrons 
which are on large distances (like 4s and Ap electrons). 
The 3d electrons are on smaller distances and have larger 
effect on atomic core. However, in most of the cases (see 
Table HI only one among six 3d electrons change its state. 
Therefore their effect on atomic core is also small. More 
detailed discussion on the effect of valence electrons on 
atomic core can be found in Refs. [5^, [3(| ■ 

All configurations in Table U correspond to an open- 
shell system. We perform the calculations stay- 
ing within central-field approximation but using frac- 
tional occupation numbers. As a result we have 
23 singe-electron basis states for valence electrons: 
3dg L, 3dg l L, 4s^\ 4p^L, 4pg l L. Here index i is the set 



number (as in Table [TJ . Note that total number of basis 
states is less than 5 times number of sets because many 
configurations don't include particular single-electron 
states. Note also that our basis set is non-orthogonal, 
e.g. < { 3 rfg 2 y 2 1 3 c?g y 2 ) < 1. The implications of this fact 
will be discussed below. 

The effective Hamiltonian for valence electrons has the 
form 



i=i i<j 



e /rij, 



(7) 



h\{ri) is the one-electron part of the Hamiltonian 



hi = ca ■ p + {[3 — l)mc 2 



Ze 2 
r 



V core + SV. (8) 



Here a and (3 are Dirac matrixes, V core is Hartree-Fock 
potential due to 18 core electrons (Is 2 . . . 3p 6 ) and 5V 
is the term which simulates the effect of the correlations 
between core and valence electrons. It is often called 
polarization potential and has the form 



SV = - 



2(r 4 + a 4 ) ' 



(9) 



Here a p is polarization of the core and a is a cut-off pa- 
rameter (we use a — as)- The form of the SV is chosen to 
coincide with the standard polarization potential on large 
distances (— a p /2r 4 ). However we use it on distances 
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where valence electrons are localized. This distances are 
not large, especially for the 3d electrons. Therefore we 
consider 5V as only rough approximation to real correla- 
tion interaction between core and valence electrons and 
treat a p as fitting parameters. The values of a p for each 
configuration of interest are presented in Table |U They 
are chosen to fit the experimental position of the config- 
urations relative to each other. The value of a p for the 
3d 6 Ap 2 configuration is taken to be the same as for the 
ground state configuration because actual position of this 
configuration in the energy spectrum is not known. For 
all configurations the values of a p are very close. This 
is not a surprise since the core is always the same. One 
can probably say that small difference in a p for different 
configurations simulates the effect of incompleteness of 
the basis and other imperfections in the calculations. 



A. CI calculations with a non-orthogonal basis. 

As it was mentioned above we have a set of single- 
electron states which is not orthogonal. The 3d, As and 
Ap states in the configurations listed in TableUare similar 
but not the same. In principle, it may lead to complica- 
tion in the CI procedure, starting from non-orthogonality 
of many-electron basis states which would lead in turn to 
complications in calculation of matrix elements and ma- 
trix diagonalization. However, most of these complica- 
tions can be avoided by appropriate selection of the con- 
figurations included in the calculations. It is sufficient to 
obey the two rules: 

• Forbid configurations which have singe-electron 
states taken from different sets, e.g. 3d™3d™4sfc4s;. 
Here i, j, k and I are set numbers as in Table|T] (i ^ j 
or/and k ^ I) and m and n are number of electrons 
in each of the 3d state (m + n = 6). 

• Don't generate additional configurations by excit- 
ing electrons to the orbitals of the same symmetry, 
e.g. 3dHs 2 -> 3d 6 4s5s. 

In present calculations we use only those configurations 
which are listed in Table U 

If all single-electron states for every many-electron ba- 
sis state are taken from the same set then the many- 
electron basis states remain orthogonal to each other. 
Indeed, states of the same configuration are orthogonal 
to each other as in the standard CI technique. States of 
different configurations are orthogonal because at least 
one electron changes its angular symmetry in the transi- 
tion between the configurations. For example all states 
of the 3d 6 4s 2 configuration are orthogonal to all states 
of the 3d 7 As configuration because of the s — d transition 
involved. 

Since many-electron basis functions remain orthogonal 
matrix diagonalization is not affected. Calculation of the 
matrix elements between states of the same configuration 
is not affected as well. The only part of the CI procedure 



which is affected is calculation of matrix elements be- 
tween basis states of different configurations. Here single 
electron part hi §5§ of the Hamiltonian does not con- 
tribute because this is a scalar operator which cannot 
change angular symmetry of a single-electron state. Only 
Coulomb integrals contribute to the matrix elements and 
their calculation must be accomplished by the product of 
overlaps between similar states from different sets. For 
example, Coulomb interaction between the 3d 6 4s 2 and 
3d 6 4p 2 configurations has the form (in non-relativistic 
notations): 

Fi(4si,4p 3 , 4s 1 ,4p 3 )(3d 1 |3d 3 ) 6 . 

Here F\ is dipole Coulomb integral, indexes I and 3 nu- 
merate basis sets as in Table HI (3di |3c?3) is the overlap 
between different 3d functions. 



III. RESULTS AND DISCUSSION 

Neutral iron is an interesting system as a challenge for 
the calculations and as a candidate for the search of the 
variation of the fine structure constant. There is strong 
configuration mixing between the 3d 6 4s 2 and the 3d 7 4s 
even configurations in the ground state and the 3d 6 4s4p 
and the 3d r Ap configurations for the odd excited states. 
The latter mixing is a fortunate feature which makes Fe I 
a convenient object for the analysis. Let us elaborate. It 
is important to have relativistic frequency shifts of the 
atomic transitions used in the analysis to be as large as 
possible. The value of the shift depends on how many 
electrons change their states in the transitions and how 
large is the change of electron momentum in each single- 
electron transition (see, e.g. [HI)- The transition between 
3d 6 4s 2 and 3d 6 4s4p configurations is basically a 4s — Ap 
transition. However, mixing with the 3d 7 4s configuration 
in the upper state adds one more single-electron transi- 
tion (4s — 3d) and makes the frequency shift larger. Note 
that the presence of both the 3d 6 4s4p and the 3d 7 Ap con- 
figurations is important. The first configuration is needed 
for the transition to the ground state to be strong elec- 
tric dipole transition, otherwise it will not be observed. 
The second configuration is needed for the relativistic 
frequency shift to be large. It is fortunate that strong 
configuration mixing between these two configurations 
takes place for most of the low odd states of Fe I. 

On the other hand this strong configuration mixing is 
a big challenge for the calculations. It makes the results 
for the relativistic energy shifts (the g-coefficients) to be 
unstable since they are very sensitive to the value of the 
mixing. Note that the configuration mixing in the ground 
state is also important. The admixture of the 3d 7 4s con- 
figuration adds the contribution of the 3d — Ap transition 
to the relativistic frequency shift. This contribution has 
an opposite sing as compared to the 4s — Ap transition 
which adds to the instability of the results. 

Since configuration mixing is very sensitive to the en- 
ergy intervals between the states the most reliable re- 



0.2 0.4 0.6 0.8 1 

(«/«o) ! 

FIG. 1: Odd-parity energy levels of Fe I with total momentum 
J = 3 as functions of the fine structure constant. 
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FIG. 2: Odd-parity energy levels of Fe I with total momentum 
J = 4 as functions of the fine structure constant. 



suits can be obtained in the calculations which reproduce 
correctly experimental spectrum. In present calculations 
this is achieved with the use of the core polarization term 
© in the Hamiltonian and fitting the data by changing 
the core polarizability parameter a p . Note however that 
only fine tuning was needed since in the end the values 
of the a p for different configurations turned to be very 
close to each other (see Table fl} . 

Another source of possible numerical instability of the 
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FIG. 3: Odd-parity energy levels of Fe I with total momentum 
J = 5 as functions of the fine structure constant. 



results for particular states is level pseudo-crossing (see, 
e.g. [H, HH). Energies of the states when considered 
as function of a 2 may come close to each other in the 
vicinity of the physical value of a. Then small error in 
the position of level crossing may lead to large error in 
the g-coefficient which is actually the slop of the curve 
E(a 2 ) (see, formula ([2]). To investigate whether this is 
the case for Fe I we plot the energies of few low odd states 
of Fe I with total momentum J = 3, 4 and 5 as function 
of a 2 from non-relativistic limit a = to the physical 
value of a. The results are presented on Figs [U [2] and 
[3] As can be seen from the pictures, there are multiple 
level crossing for states with J = 3 and J — A. However, 
all these crossings take place on safe distance from the 
physical value of a and are very unlikely to cause the 
instability of the results. Another interesting thing to 
note is that the energies are practically linear functions 
of a 2 in all cases. 

Table HI] presents experimental and theoretical energies 
and (/-factors of the lowest even states of Fe I. The g- 
factors are useful for the identification of the states and 
for control of configuration mixing [l6l |. As can be seen 
the experimental data are reproduced in the calculations 
with very good accuracy for both the 3<i 6 4s 2 and 3d 7 4s 
configurations. 

Table IIIII presents experimental and theoretical ener- 
gies and g-factors of the lowest odd states of Fe I in which 
the 3d 6 As4p configuration dominates. Theoretical rela- 
tivistic frequency shifts (q-coefficients) are also presented. 
The q-coefficients were obtained by numerical differentia- 
tion using formula (j4|). Note that only states with J=3,4 
and 5, for which electric dipole transition to the ground 
state is possible are needed for the analysis. However, we 
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TABLE II: Energy levels (cm 1 ) and g- factors of the lowest 
even states of Fe 



Conf. 



Term J 



Experiment" 



Calculations 









Energy 


q 


Energy 


a 


3d" 4s' 1 a 


b D 


4 


0.000 


1.500 





1.4995 






3 


415.932 


1.500 


464 


1.4997 






2 


704.004 


1.500 


790 


1.4998 






1 


888.129 


1.500 


1000 


1.4998 









978.072 




1103 


0.0000 


3d 7 4s a 


5 F 


5 


6928.266 


1.400 


6862 


1.3996 






4 


7376.760 


1.350 


7374 


1.3496 






3 


7728.056 


1.249 


7779 


1.2497 






2 


7985.780 


0.999 


8078 


1.0000 






1 


8154.710 


-0.014 


8275 


0.0010 


3d 7 4s a 


3 F 


4 


11976.234 


1.254 


13040 


1.2496 






3 


12560.930 


1.086 


13702 


1.0835 






2 


12968.549 


0.670 


14171 


0.6676 
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TABLE III: Energy levels (cm 1 ), p-factors and relativistic 



energy shifts (g-factors, cm 
configuration of Fe. 



for the states of of the 3d 4s4p 



Term J 



'D° 



7po 



7po 



Experiment" 



Calculations 





nil nrriT 


9 




9 


9 







i ^Q7 


±y 100 


1 ^Q87 


/1QO 


A 

4 


lyODZ.44U 




1 QQQA 


1 p. A Qfl 


fifi9 
DDZ 


3 


19757.033 


1.746 


19611 


1.7485 


891 


2 


19912.494 


2.008 


19793 


1.9976 


1092 


i 
i 


9nni q 
zuuiy .doo 


9 QQQ 


1 QQ91 


z.yyou 


1 937 


/• 




OOCCH AO! 

22DO0.421 


i /i no 


o i aao 
zlDo3 


i.4yy ^ 


K CO 

082 


5 


22845.868 


1.498 


21891 


1.5005 


827 


4 


22996.676 


1.493 


22062 


1.5026 


982 


3 


23110.937 


1.513 


22189 


1.5029 


1103 


2 


23192.497 


1.504 


22282 


1.5026 


1184 


1 


23244.834 


1.549 


22338 


1.5029 


1227 





23270.374 




22366 


0.0000 


1246 


4 


23711.457 


1.747 


22543 


1.7470 


491 


3 


24180.864 


1.908 


23034 


1.9136 


983 


2 


24506.919 


2.333 


23440 


2.3309 


1316 



present ^-coefficients for all states for better illustration 
of the accuracy of the calculations. In the linear in a 2 
approximation the difference in g-coefficicnts for states 
of the same fine-structure multiplet is equal to the fine 
structure interval between this states. As can be seen 
from Figs. 1 1121 and [3] the dependence of the energies on 
a 2 is very close to linear indeed. Therefore, comparing 
the data for the fine structure and q is another test of 
the calculations. 

Table HVl presents the data similar to those of Table ILTTl 
but for the states where the 3d 7 4p and 3d 5 4s 2 4p configu- 
rations dominate. The values of the ^-coefficients for the 
states of the 3d 7 4p configuration are larger than those 
of the 3d 6 4s4p configuration. This is due to additional 
contribution from the 4s — 3c? singe-electron transition as 
it was explained above. 

It is interesting that similar to the case of the ion 
Fe II [III neutral iron also has some negative shifters (q < 
0). Corresponding states belong to the 3d 5 4s 2 4p con- 
figuration. Negative sign of q is due to the dominant 
contribution from the 4p — 3d single-electron transition. 
The data are presented in Table HVl Note however that 
the spin of these states is different from the spin in the 
ground state. This means that the electric dipole transi- 
tion is suppressed by conservation of spin and goes only 
due to relativistic effects. This in turn probably means 
that the transitions may be too weak to be observed. 

We estimate the accuracy of present calculations of 
the ^-coefficients to be on the level of 20 to 30%. The re- 
sults were obtained with a very simple method which uses 
small number of basis functions and some semi-empirical 
fitting. The main challenges for more accurate calcula- 
tions are strong configuration mixing and large number 
of valence electrons. Further development of the methods 
or the use of supercomputers might be needed for better 
accuracy of the calculations. 



z 


5 D° 


4 


25899.987 6 


1.502 


26428 


1.4979 


999 






3 


26140.177 


1.500 


26679 


1.4984 


1223 






2 


26339.691 


1.503 


26924 


1.4976 


1450 






1 


26479.376 


1.495 


27094 


1.4971 


1616 









26550.476 




27174 


0.0000 


1705 


z 


5po 


5 


26874.549 6 


1.399 


27432 


1.3999 


880 






4 


27166.819 


1.355 


27702 


1.3517 


1180 






3 


27394.688 


1.250 


27947 


1.2530 


1402 






2 


27559.581 


1.004 


28119 


1.0041 


1568 






1 


27666.346 


-0.012 


28213 


0.0062 


1680 


z 


5po 


3 


29056.321 6 


1.657 


29340 


1.6643 


859 






2 


29469.020 


1.835 


29795 


1.8307 


1310 






1 


29732.733 


2.487 


30118 


2.4966 


1594 


z 


3po 


4 


31307.243 


1.250 


32356 


1.2504 


1267 






3 


31805.067 


1.086 


32883 


1.0885 


1808 






2 


32133.986 


0.682 


33263 


0.6767 


2177 


z 




3 


31322.611 


1.321 


32032 


1.3314 


1456 






2 


31686.346 


1.168 


32464 


1.1662 


1843 






1 


31937.316 


0.513 


32750 


0.5035 


2119 
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IV. CONCLUSION 

We have calculated relativistic frequency shifts for a 
number of the lower odd states of Fe I. Some of these 
states were observed in the quasar absorption spectra. 
Calculations show that due to strong configuration mix- 
ing the values of the shifts are large and vary significantly 
between the states. This makes Fe I to be a good can- 
didate for the search of variation of the fine structure 
constant in quasar absorption spectra. 
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TABLE IV: Energy levels (cm - ), g-factors and relativistic 
energy shifts (q-factors, cm -1 ) for the states of the 3d 7 4p and 
3d 5 4s 2 4p configurations of Fe 



Conf. 


Term 


J 


Experiment" 


Calculations 








Energy g 


Energy 


9 


q 


3d Y 4p 


y b D° 


4 


33095.937" 1.496 


32680 


1.4511 


2494 






3 


33507.120" 1.492 


33134 


1.3492 


3019 






2 


33801.567 1.495 


33466 


1.1053 


3423 






1 


34017.098 1.492 


33705 


0.1658 


3754 









34121.58 


34007 


0.0000 


3723 



3d 7 4p y 5 F° 



5 


33695.394* 


1.417 


32522 


1.3964 


2672 


4 


34039.513 


1.344 


33029 


1.3913 


3021 


3 


34328.749 


1.244 


33404 


1.3881 


3317 


2 


34547.206 


0.998 


33705 


1.3777 


3536 


1 


34692.144 


-0.016 


33909 


1.3347 


3678 



3d 7 4p z 5 G° 5 34782.416 

6 34843.94 

4 35257.319 

3 35611.619 

2 35856.400 



1.218 33978 

1.332 33687 

1.103 34363 

0.887 34661 

0.335 34883 



1.2487 3024 
1.3330 

1.1406 3520 

0.9153 3864 

0.3505 3464 



3d 7 4p z S G° 5 35379.206 1.248 34506 1.2209 3340 
4 35767.561 1.100 35042 1.0731 3697 
3 36079.366 0.791 35474 0.7671 4096 



3d 7 4p 



3d 5 4s 2 4p y 7 P° 



4 


36686.164 


1.246 


35697 


1.2425 


3085 


3 


37162.740 


1.086 


36227 


1.0863 


3487 


2 


37521.157 


0.688 








2 


40052.030 


2.340 


40529 


2.3278 




3 


40207.086 


1.908 


40677 


1.8883 


-2472 


4 


40421.85 


1.75 


40926 


1.7491 


-2287 



a NIST, Ref. HI 

"States observed in quasar absorption spectra 
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